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Sagdeyev and Shafranov have shown that the absorption of microwaves in a hot 
plasma in a steady magnetic field can be calculated in simple closed form with the help of 
the Boltzmann equation, provided that the effect of collision can be ignored. 

The present paper is restricted to the special case of propagation of circularly polarized 
waves parallel to the magnetic field, and the extraordinary ray, in magneto-ionic termi- 
nology, is given special attention. It is shown that the formula given by Sagdeyev and 
Shafranov for this case can be deduced by considering the motions of individual electrons 
by elementary dynamical methods, using the concepts of Doppler shift and velocity distri- 
bution functions to obtain a macroscopic conductivity formula for a high-temperature 
plasma. From this, the absorption is easily calculated. 

It is emphasized that the calculation in no way depends upon the assumption of a 
Maxwellian velocity distribution function. The absorption can in fact be obtained in 
closed form for any arbitrary velocity distribution function. 

This suggests that a diagnostic technique for the determination of velocity distribution 
could be based on measurements of absorption of the extraordinary ray, and the potential- 
ities and limitations of this proposal are briefly discussed. 



1. Introduction 

Sagdeyev and Shafranov 2 have shown that the 
absorption of microwave energy in a hot plasma can 
be calculated in simple closed form from the Boltz- 
mann equation provided that the effect of collisions 
can be ignored. 

In the present paper, the same formula is obtained 
from a simple dynamical-collision al argument, and 
it is shown that if the collision frequency is small in 
comparison with the Doppler frequency shift due to 
thermal motions of the electrons, the value of the 
collision frequency has only a second-order effect on 
the absorption coefficient. Furthermore, the results 
of Sagdeyev and Shafranov are extended to cover 
the case of an arbitrary, rather than a Maxwellian, 
velocity distribution function, and a microwave 
diagnostic technique for determining the distribution 
functions in low density hot plasmas is proposed. 

2. Equations of Motion 

In these equations, the effect of the radiofrequency 
magnetic field will be ignored. The electric field is 
assumed in the first instance to have the form 



E X =E sin o)t 
E y =E cos cot 



} 



a) 



i Contribution from Microwave Laboratory, W. W. Hanson Laboratories of 
Physics, Stanford University, Stanford, Calif., while author was on leave from 
the Department of Electrical Engineering, The University of Sheffield, Sheffield, 
England. The research reported in this paper was sponsored by the Air Force 
Cambridge Research Center, Air Research and Development Command, 
Bedford, Mass., 1959. 

- Sagdeyev and Shafranov, Absorption of high-frequency electromagnetic 
field energy in the high-temperature plasma. Proceedings of the Second Inter- 
national Conference on Peaceful Uses of Atomic Energy, Geneva, September 
1958. Paper P/2215. 



It is also assumed that there is a constant mag- 
netic field along the z-axis of value B {) . If the 
strength of this field is expressed in terms of the 
cyclotron frequency for electrons u c =eB lm, the 
equations of motion can be written 



sin Cot^WcVy 



dVx = ___e_E^ 
dt m 

dv y eE , 

-77= cos <at±w c v x 

dt m 



(2) 



The upper sign corresponds to a magnetic field di- 
rected along the positive 2-axis. Looking in this 
direction, the electric vector described by (1) rotates 
anticlockwise. The natural direction of rotation of 
electrons in this magnetic field is clockwise. Thus 
if (1) described the electric field of a circularly polar- 
ized wave traveling along the 2-axis in the positive 
direction, the upper sign will correspond to the 
ordinary ray, and the lower sign to the extraordinary 
ray, in ionosphere theory terminology. 

The extraordinary ray is the case of interest here, 
for this wave can be heavily absorbed by synchronous 
acceleration of the electrons by a cyclotron-type 
mechanism. In what follows, the lower sign will 
be taken. To use the equations for the ordinary 
ray, it is only necessary to change the sign of u c . 

From eq (2) we find 






dt' 

d\ 
dt 2 



m 



eE, 



+<a 2 c Vy=+— ^ (co+co c ) sin oot 



m 



(3) 
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3. Collisions and Boundary Conditions 

For our present purpose the details of the collision 
process are not important, since it will emerge that 
the collision frequency does not enter into the final 
formula for absorption. 

We shall assume that an electron is brought to 
rest at each collision. Thus, for an electron which 
makes a collision at time t', its motion until the next 
collision takes place can be found by solving the 
eq (3) subject to the conditions 



~-Vy=0 



dv, 
~dt 



-rf= = -sinco^ 



m 



(4) 



dt 



eE ., 

= cos cot 

m 



The last two formulas follow by substituting the 
first in (2). The solution of (3) appropriate to the 
boundary conditions (4) is given by 



eE d 



Vy = 



m(co— co c ) 
eE 



m(co — w c ) 



[cosco^ — cos{co c Z+(w— co c )t'}] (5) 

[sincot— sin {co c t-\-(co — co c )t f }]. (6) 



Equations (5) and (6) can be expressed more com- 
pactly thus 



%—3V v = 



eE 



m(co — co c ) 



\ e i ut e H « c H-(«- 



*n (7) 



Equation (7) gives the x and y components of 
velocity at time t of an electron which was at rest 
at time V . 

Assuming that n, the number of electrons per 
unit volume, is sufficiently large, we can assume 
that, on the average, at any time t the number of 
electrons per unit volume whose most recent collision 
occurred between t' and t fJ rdt' is given by 



dn—nve~ 



v(t-t') 



dt'. 



(8) 



Using this formula, which in fact defines the collision 
frequency v, we can calculate the current density 
at any time t. The contribution to J x from the 
electrons whose most recent collision was in the 
interval t' to t fJ rdt f is dJ x =—nev x dn, where v x is 
given by (5). Hence, using (7), we find 



Jx—jJy 



—nve 2 E 



m(a 



■.) 



f 

t/ — 



p3ut_ 



o 0\ ut+(u 



")]e- 



v(t-t') 



dt'. (9) 



Carrying out the integrations in (9) and separating 
the real and imaginary parts, we find 



j x =- 



ne 2 E Vv sin ut— (co— co c ) cos cot 



m 



■[ 



» 2 +(<o- 



cY 



] 



(10) 



Jy : 



ne 2 E [~v cos cot-\-(co— co c ) sin cot 



m 



I s 



} 



4. Effect of Velocity Distribution 



(11) 



So far, the formulas we have derived are familiar 
in magneto-ionic theory. We now consider the effect 
of thermal velocities of the electrons. If the electric 
field we have been discussing is associated with a 
wave of frequency co traveling in the direction of 
increasing z, with phase constant /3, then an electron 
traveling in the same direction with velocity v z will 
be subjected to a field of angular frequency 
co — $Vz=co say. 

If dn electrons/unit volume have a ^-component 
of velocity between v z and v z -\-dv z , and if the velocity 
distribution is Maxwellian, we have 



dn =n 



/ m V' 
WkTj 



-mv\l2kT , 



(12) 



Equations (10) and (11) can be used to calculate 
the contribution to the current density components 
J x and J y in a frame of reference traveling with the 
electrons, i.e., with velocity v z along the original 
2-axis, if we substitute dn for n, and interpret co as 
ooo—fiVz- Thus 



dJ x = 



n e 2 E Q 



?n 



/ m V 



-mv 2 l2kT 



sin cot— (oo— co c ) cos cot 



,*+ («-«<)' 



] 



dv t - (13) 



Transforming this contribution to the current back 
to the original frame of reference does not affect its 
amplitude but restores the frequency to the original 
value co - Thus, summing all such contributions 



J x = 



n e 2 E 



-mv 2 J2kT 



m 



Similarly, 

n o e 2 E 



( m Y /2 f H 

tv sin co t — (co — co c ) cos co t"l 



Jy= 



-mv 2 J2kT 



dv z - (14) 



m 



/ m Y 72 f + 
WkTj J_ o 

v cos co t-\-(co— co c ) sin co t 



D 



^+(a 



-co,) 2 



] 



dvs (15) 



Note that parts of J x and J y which are in phase 
with the electric field components E x and E y are 
equal, and correspond to a transverse conductivity 
ci given by 



n e 2 / m \ 1/2 f + 



-mv 2 /2kT 



h 



+ ("o- 



•] 



dv,. (16) 
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Similarly, the transverse permittivity e Y is 



«j = 



mv z /2kT 



n e 2 / m \ 1/2 f + 
" t0 i^fflW J_„ 

r ("o— w, 



-flO 



v 2 +(«o— We— j8»«) 



»]<*>.. (17) 



The effect of collisions on the t? g component of ve- 
locity has not been considered, since the absorption 
phenomenon is not affected by this process. We do 
not need to have detailed knowledge of v z for each 
electron, we only need to know the number of elec- 
trons in any given small range of velocities. 

5. Wave Propagation in a Hot Plasma 

We assume now that the plasma temperature T 
is large 3 and the collision frequency small, so that 



p<0- 






(18) 



In this case the bracketed factor in the integrand of 
(16) is small unless v z is such that 



w — w c — fto z <Cy. 



(19) 



Within this small range of values of v z the ex- 
ponential factor can be given the value corresponding 
to v z =(co — co c )//3, and taken outside the integral 
sign. It is then a simple matter to evaluate the 
integral, and the resulting formula for the conduc- 
tivity is 



_ Trn () e 2 / m \ 1/2 



-m/2kT[( UQ - Ue )/p] 2 



(20) 



If the conductivity is small, the absorption coeffi- 
cient k can be calculated very simply from the 
formula 



a cr± 



*± 



2a*o2> 2 



(21) 



where a and /3 are the attenuation and phase con- 
stants respectively, and p is the real part of the 
refractive index. If, in accordance with Sagdeyev 
and Shafranov, we define V\\ by the equation 



we find 



Vtt 



C (Ot. 



K -2^|,o,^ eXP 



V m 

\ v\\Y\ 



OJ 



)'}■ 



(22) 



(23) 



| Here co p is the plasma frequency (oo p 2 =ne 2 /meo), and 
co and co c are the frequency of the wave and the 
electron cyclotron frequency respectively. 

This formula agrees exactly with the result of 
Sagdeyev and Shafranov, as stated in eq (13) of 
I their paper (see footnote 2). 

The calculation of ej_ from (17) is carried out in 



3 Since this is the only temperature involved, we use T for Til, the effective 
temperature associated with thermal motion of electrons parallel to the magnetic 
field lines. 



the appendix. When this is done, we find 



V\\ U%P 



^[(uo-coc)^] 2 f c / q? — coA 

-&(*?)■+ -}■ <m) 



This also agrees with the form given by Sagdeyev 
and Shafranov, although their result is expressed in 
terms of an error integral. (The lower limit of this 
integral, which is missing, is presumably zero.) 
There is also an error in the sign of the second term 
in their formula. 

6. A Microwave Diagnostic Technique for 
Determining the Electron Velocity Distri- 
bution in a Hot Plasma 

It is easy to see that the method employed for 
integrating (16) does not depend in any way on the 
velocity distribution. This need not be Maxwellian, 
but can be quite arbitrary. 

Let the velocity distribution parallel to the mag- 
netic field be defined by formulas 

J(v z )dv z =l (25) 

in place of the Maxwellian distribution formula eq 
(12). 

We find for the conductivity of the plasma the 
very simple expression 

iMoe 2 ,/co — coA 
The attenuation coefficient can be written 



7T6> C0 



W<) rA>0 — «A 



2m & 



Thus the distribution function can be expressed in 
terms of experimentally observable quantities as 
follows 

^(^y^«jt. (26) 



7T6 2 /Xo 



w 



If the right-hand side of (26) is plotted as a function 
of (co — o) c )lPj the electron velocity distribution par- 
allel to the magnetic field is obtained. The area 
under this curve will give the electron density. In 
practical applications, the distribution function may 
depend on z. If this is the case, the total attenuation 

measured is approximately adz along the transmis- 
sion path, provided that the rate of change of 
attenuation per wavelength of path is not too great. 
Thus from the total attenuation we can determine, 
using (26), the average distribution function, the 
average being taken along the microwave path. 
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7. Wave Propagation in a Cold Plasma 

For mulas for a± and ej_ in a cold plasma, in which 
v^>/3^kT/m, are easily obtained from (16) and (17) 
if the collision frequency can be assumed independent 
of v z . The familiar magneto-ionic theory formulas 
are obtained: 



0-_L=- 



ne' 



m ^ 2 +(w -co c ) 2 



e 



(cop— 0) c ) 
e m v 2 +(u — co c )2 



>• 



(27) 



Note that if ^ ^>(co — co c ), o-_l is proportional to 1/v, 
whilst if j><C(co — co c ), cjL is proportional to p. 

8. Discussion 

The purpose of this section is to contrast the be- 
havior of the conductivity of hot and cold plasmas 
with respect to variation of collision frequency. 

Consider first the cold plasma conductivity as 
described by (27). For a large collision frequency, 
the conductivity is proportional to \\v. The physical 
meaning is clear. It is easy to show from (7) that 
an electron starting from rest at time t~0 acquires 
in time t a kinetic energy given by 



W= 



le 2 E 



-[(?)~ 



m 



(28) 



assuming that no collisions take place. Differenti- 
ating, we find 



dW__e 2 E I sin (co-co c )£ 



dt 



(co— co c ) 



(28a) 



If the collision frequency v is large, in comparison 
with (o) — o) c ), the time interval t, during which the 
gain in kinetic energy described by (28a) can take 
place before the process is interrupted by another 
collision, is small, and (co — co c ) £<Cl. Hence, (28a) 
can be replaced by 



dW^e 2 El 
dt m 



(29) 



The average rate of absorption per electron is given 

by 



dt m v 



(30) 



Thus the physical significance of the inverse varia- 
tion of conductivity with collision frequency, when 
the collision frequency is large in comparison with 
the frequency difference, co — co c becomes clear. 



When the collision frequency is small in comparison 
with (co — w c ), the amplitude of oscillation, and hence 
the mean kinetic energy acquired by the electrons, 
is not limited by collisions, but by the lack of syn- 
chronism between the frequency of the alternating 
force by the wave and the natural frequency of 
gyration of the electrons. If co >co c , the amplitude 
of oscillation is limited mainly by electron inertia, 
while if co <co c , it is limited by the obstructing effect 
of the magnetic field. If (28) is averaged over the 
long time which elapsed between collision, we find 

Try e*El 



W-. 



"m(co - 



•cY 



If it is assumed that, on the average, this energy, 
or a definite fraction of it, is given up by an electron 
at each collision, the average rate of absorption of 
energy per electron is given by 



vW= 



e 2 El 



m(co — co c ) 2 



(31) 



and is clearly directly proportional to v, the collision 
frequency. Thus the reason for the variation of 
conductivity with collision frequency predicted by 
(27) is clear, and the physical mechanism by which 
the energy is absorbed has been exhibited. Energy 
extracted from the wave appears first as kinetic 
energy of the electron, and is then transferred to 
molecules or ions of the gas by collisions. This 
absorption process is called collisional absorption, and 
as we have seen, depends strongly on the collision 
frequency. 

In direct contrast to this situation, the absorption 
by a "hot" plasma, in which ix^Q^kT/m, is inde- 
pendent of v to first order, as eq (20) shows. The 
reason for this is that in a hot plasma, the effective 
frequency seen by an electron depends on its veloc- 
ity in the direction of the propagation of the wave, 
due to the Doppler shift. Since the electron veloci- 
ties vary widely in a hot plasma, a wide range of 
effective frequencies exists, and the hot conductivity 
is essentially an average of the cold conductivity 
over all frequencies. That this is independent of v 
can be easily seen using (27) 



f + 00 . ne 2 f J 
J-co mj- 



2 +("-co c ) 



do. 



m^C^ v_ 

?nj -oo V 2 + < 



, dw x = 



Tne" 
m 



(32) 



Provided that v is small enough, the integrand is 
very small except when (co— w c )~*>, so that a very 
small range of effective frequencies (and electron 
velocities) contributes significantly to the absorption. 
Provided that the plasma temperature is high 
enough, the variation of the distribution function 
over this small range of significant electron velocities 
may be neglected, and the distribution function can 
be regarded as constant as was done in obtaining 
eq (20). 
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Further, if v depends on r z , the analysis is still 
valid, provided that the dependence is not strong 
in the small range of significant velocities. 

The case of no collisions at all can be regarded as 
(he limiting ease v— >0; the "cold conductivity" 
term in the integrand of (16) then has the character 
of a delta function, and eq (10) is exact. 

Thus, the agreement of our analysis with that of 
Sagdeyev and Shafranov is to be expected, in spite 
of the fact that in their treatment it was assumed ab 
initio that no collisions occurred, whereas the 
present analysis necessarily involves the considera- 
tion of collisions. 

9. Conclusions 

Summarizing, the present analysis has shown: 

(a) That in the absorption of microwaves by a hot 
plasma in a magnetic field, collisions have a second- 
order effect only; 

(b) That the results of Sagdeyev and Shafranov, 
obtained by solving the Bolt zm aim equations for 
the distribution functions by a perturbation technique, 
neglecting collisions at the outset, can be obtained by 
a direct ballistic analysis of the motion of individual 
electrons ; 

(c) That the assumption of a Maxwellian velocity 
distribution function made by Sagdeyev and 
Shafranov is not necessary to the analysis, and that 
a closed-form solution for conductivity can be 
obtained for any arbitrary distribution function : 

(d) That a study of the variation with frequency of 
the attenuation and phase constant of a circularly 
polarized electromagnetic wave in an ionized plasma 
in a magnetic field can in principle yield the velocity 
distribution function for the thermal motion of 
electrons parallel to the magnetic field lines. There 
may be considerable difficulty in applying this 
technique in some cases. For example, the extra- 
ordinary ray only must be employed, and a high 
degree of discrimination against the ordinary ray in 
the launching and/or receiving antennas will be 
necessary if errors due to the ordinary ray, which is 
attenuated much less strongly, are to be negligible. 
Also diffraction effects will be serious except at 
microwave frequencies, so that the method can only 
be applied in the simple form suggested here if the 
magnetic field strength is of the order of 3,000 gauss 
or more, so that the cyclotron frequency falls in the 
microwave-frequency band. 

As an alternative to attempting to eliminate the 
ordinary ray entirely, one could start with a linearly 
polarized wave and measure the ellipticity of the 

\ received signal. This would have some practical 
advantages, and the necessary theory could be 
developed very simply from the equations given here. 
A more serious limitation arises when the order of 
magnitude of the attenuation at resonance is esti- 
mated. For n= 10 11 electrons/cm 3 T=10 6 °K, and 
fo= 10,000 Mc/s, the attenuation coefficient at 

} cyclotron resonance has a value of about 100 db/cm 
or the attenuation length is about 1 mm. For 
densities of 10 8 electrons/cm 3 under the same con- 
ditions, however, the attenuation coefficient is about 



0.1 db/cm, and for such a low density plasma the 
method might be feasible. 

It should be noted that the effect of the radio- 
frequency magnetic field has not been considered in 
the present analysis. While this seems unlikely to 
lead to serious errors at electron temperatures of a 
million degrees or so, these effects might be ap- 
preciable for temperatures of the order of a hundred 
million degrees, when the electron velocity is only 
one order of magnitude less than the velocity of 
light. 

No comparison has been made with the more 
elaborate theory of Drummond 4 , which starts from 
the Boltzmann equation but retains the collision 
term, for which a suit a hie approximation is later 
introduced. It would be very desirable to make 
such a comparison, which would probably shed 
further light on the significance of collisions in a 
practical situation. 

10. Appendix 
We have to evaluate 

/= f + " e -« !/*■ fo-u-fto.) d (1) 



Let t = jm/2kTv z : t =-\ l m/2kT[(w — <o c )/#|. Let us 
also put v=0. We get 



C +m e-' 2 dt 
to-t' 



Hi 



Now put z=t—t . 



1 C +m e" (2+ 'o )2 



/=— dz. 



(2) 



(3) 



Because of the singularity at 2=0, we interpret / in 
the following way 

-jSJ^-fllim-f f"Y« , -»o-^+ r + V« f -*««^y 

Replacing z by — z in the first integral, and com- 
bining the two integrals gives 



+- 



Jo 



-sinh 2t z dz. 



(4) 



This integral is easily evaluated by expanding 
sinh 2t z and integrating term by term. This leads to 
the result given in eq (24) . 
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4 J. E. Drummond, Basic microwave properties of hot magnetoplasma, Phys. 
Rev. 110 (April 15, 1958). 

(Paper 64D5-89) 
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